function J = computeCost(x, y, theta)
    m = length(y);
    % J = 0;
    % J = sum((x * theta - y).^2) / (2 * m);
    h = x * theta - y;
    J = h' * h / (2 * m);

end
